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Emission line variability of RS Ophiuchi * 



in 
o 
o 

(N 



> 

m 
m 
\o 
o 
in 
o 

6 

CO 
C3 



R. K. Zamanov 1 ' 2 ! M. F. Bode 1 , N. A. Tomov 3 , J.M. Porter 1 

1 Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Birkenhead, CH41 1LD, 

2 Institute of Astronomy, Bulgarian Academy of Sciences, 72 Tsarighradsko Shousse Blvd., 1784 Sofia, Bulgaria 

3 Institute of Astronomy, National Astronomical Observatory Rozhen, POBox 136, 4700 Smolyan, Bulgaria 



Accepted 12 July 2005. Received 2005 June 6 



UK 



ABSTRACT 

We report that the Ha emission line of RS Oph was strongly variable during our 2004 
observations on a time scale of 1 month. The line consisted of both a double peaked 
central narrow component (FWHM~220 km s" 1 ) and a strongly variable broad one 
(FWHM>2000 km s" 1 ). The base of the Ha line was very broad with FWZI«4600 
km s _1 on all spectra from 1986 to 2004. The variability of the broad component 
extends from —2000 to +2000 km s _1 . Most probably this is due to either blobs 
ejected from the white dwarf (with a typical blob mass estimated to be -10" 10 M ) 
or a variable accretion disk wind. We also detected variability of the HeIIA4686 line 
on a time scale shorter than 1 day. The possible origin is discussed. 
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1 INTRODUCTION 

RS Ophiuchi (HD 162214) is a recurrent nova which under- 
went its last major outburst in 1985 (see Bode 1987 and 
references therein). Following Dobrzycka & Kenyon (1994), 
RS Oph has a binary period of 460 days with orbital in- 
clination to the line of sight of 30° — 40°, and spectro- 
scopic ephemeris T = JD2444999.9 ± 29.3 + 460 ± WE. 
The mass donor is a red giant of spectral type K4-M0 
(Miirset & Schmid 1999). The nature of the hot component 
is unclear. Most probably it is a massive white dwarf with 
Mwd -1.4 M , accreting at M >10~ 8 M Q yr" 1 (Hachisu 
& Kato, 2000). Alternatively, it has been proposed to be 
a B-type shell star with highly variable luminosity, which 
occasionally displays blue-shifted absorption features (Do- 
brzycka et al. 1996) . 

Broad, variable, emission components of the hydrogen 
lines were detected by Iijima et al. (1994), Anupama & 
Mikolajewska (1999), and Tomov (2003). Variability of the 
He lines on a time scale of hours was detected by Sokoloski 
(2003). The origin of these changes is a mystery. Here we 
report that similar variability is also evident in our 2004 
observations and discuss its possible origin. 



* based on data from ESO (program 073.D-0724) and the ING 
Archive 

f e-mail: rkz@astro.bas.bg; mfb@astro.livjm.ac.uk; to- 
mov@astro.bas.bg; 



2 OBSERVATIONS 

We secured 8 spectra in 2004 using the ESO, La Silla, 
2.2m telescope and the FEROS spectrograph. FEROS is 
a fibre-fed echelle spectrograph, providing a resolution of 
A/ A A =48000, wide wavelength coverage from about 4000 A 
to 8000 A in one exposure and a high detector efficiency 
(Kaufer et al. 1999). 

Additionally, we retrieved 4 spectra covering the region 
of the Ha line from the archive of the Isaac Newton Group 
of telescopes (ING) . A journal of observations and the mea- 
sured quantities are given in Table 1. 



3 SPECTRAL VARIABILITY 
3.1 Ha emission line 

The variability of the Ha line is immediately apparent 
(see Fig.l and Table 1). There are changes in the blue 
and red peaks, as well as in the wings of the line. On all 
spectra the base of the Ha line is very wide, the mean 
full width at zero intensity FWZI = 4650 ± 300 km s" 1 
and the mean equivalent width of the Ha emission line 
W= 112±11 A, which is about lxlO" 11 erg cm" 2 s" 1 . The 
ratio of the intensity of the blue and red peaks (V/R) varies 
from 0.51 to 0.87. Our data suggest that the Ha profile of 
RS Oph is purely in emission and consists of both a double 
peaked central narrow component (W~90 A, FWHM~200- 
250 km s _1 ) and a strongly variable broad one (W~20 A, 
FWHM>2000 km s" 1 ). 
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Table 1. Ha and HcIIA4686 observations of RS Oph. In the table are given number, date of observation, Julian Day, orbital phase, 
radial velocity of the blue peak, central dip and red peak, full width at zero intensity (FWZI) of Ha, equivalent width of Ha, equivalent 
width of HcIIA4686 , the ratio between intensities of the blue and red peak, the radial velocity of the additional emission (V emm ), and 
the origin of the spectrum. The typical errors are ±1 km s- 1 for V Wue , V dip and V red , ±250 km s" 1 for FWZI, ±5% for W Ha , ±20% 
for W 4 686, ±3% for the V/R ratio, and ±150 km s" 1 for V emm- The last column indicates the origin of the observations. 
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Figure 1. The Ha and the HcIIA4686 emission lines of RS Oph (FEROS observations). The number against each spectrum refers to the 
epoch of observation as given in Table 1. All spectra are normalized to the local continuum. 

a) The Ha central double peaked component with FWHM ~5A (220 km s~ 1 ). 

b) The same spectra as in a) but on different scales, so the Ha broad component is visible. 

c) The HeIIA4686 line. The spectrum N:8 was obtained 10 minutes after N:7, and N:9 - 24 hours later. These data demonstrate that the 
line is variable on a time scale ^1 day. 
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Figure 2. a) The RS Oph Ha profiles (normalized to the local 
continuum) for 7 different nights. The numbers again correspond 
to those in Table 1. b) After subtraction of the minimum spec- 
trum, the additional emissions are visible as bumps with 7 
A. The arrows indicate their measured radial velocities as given 
in Table 1. 



When the spectra are obtained in the same night or 
in consecutive nights, the Ha profiles are effectively identi- 
cal. However, changes are visible when the time difference is 
~1 month. To explore in more detail the velocity structure 
of the variability of the wings, assuming that the variabil- 
ity arises purely from changes in emission, we subtracted at 
each epoch a minimum spectrum (artificially generated by 
taking the minimum value from all normalized spectra at a 
given wavelength). The residuals are plotted in Fig. 2b. As 
can be seen, the variability of the wings extends from —2000 
to +2000 km s~\ with the most prominent additional emis- 
sion around —1200 and +1200 km s _1 . On each spectrum 
we measured the radial velocity of the additional emission, 
which has equivalent width W^7 A. The data are given in 
Table 1 along with other parameters. 

A variable broad component was not detected in RS 
Oph's sister system T CrB (see Stanishev et al. 2004; Za- 
manov et al. 2005). A similar component appeared how- 
ever in the 1994 observations of CH Cyg (Tomov et al., 
1996) and the velocity range of ±2000 km s _1 , was close 
to that of RS Oph. However, the time scale of appear- 
ance/disappearance was considerably shorter (~2 hours). 



3.2 Time scale of HeIIA4686 variability 

Our spectra with numbers 7, 8 and 9 show that the 
HeIIA4686 line is also variable (Fig.lc and Table 1). The 
spectrum N:8 was obtained 10 minutes after N:7, and N:9 - 
24 hours later. This suggests that the intensity of HeIIA4686 
changes on a time scale of minutes-hours (note that we did 
not detect Ha variability on such a short time scale). 

The HeIIA4686 line indicates high excitation conditions 



in the gas. In the case of RS Oph (see the next section) 
this can be related to a region located close to the bound- 
ary layer between the accretion disk and the white dwarf 
(Sokoloski 2003). The changes of the intensity of this line 
indicate changes in the region of its appearance which, in 
turn, is determined most probably by variable luminosity 
or/and temperature of the boundary layer. 



4 DISCUSSION 

Our data suggest that the Ha profile of RS Oph is purely in 
emission and consists of both a double peaked central narrow 
component and a strongly variable broad one. The Balmer 
emission lines, appearing in the extended envelopes of sym- 
biotic stars usually have an ordinary nebular profile with 
typical FWHM~100-150 km s~\ for example AG Dra (To- 
mova & Tomov 1999). FWHM increases to 200 km s _1 only 
during the active phases. The basic mechanism determin- 
ing their width is turbulence in the gas. The central narrow 
emission of RS Oph is very similar to the double peaked Ha 
line of T CrB and the two stars have practically the same 
FWHM of about 200-250 km s" 1 . The Ha line of T CrB 
is supposed to be formed in the outer part of an accretion 
disk around the hot compact object (Stanishev et al. 2004). 
By analogy with this star we could assume that the narrow 
Ha component of RS Oph is formed also around its compact 
object. This assumption, however meets some difficulty with 
the fact that the peak of this component exceeds the level of 
the continuum by a factor of 20-25 like purely nebular lines 
of the symbiotic stars. This peak height differs from that 
one of T CrB, which exceeds the level of the continuum by 
a factor of 2-6 only. 

While the appearance and the variability of the nar- 
row component is relatively common in symbiotic stars (see 
also Ikeda & Tamura 2004 and the references therein), ap- 
pearance and variability of broad component is detected in 
about 10 objects up to now. Fast (~1000 km s _1 ) bipolar 
winds/jets were detected in Hen 3-1341 and StHa 190 (To- 
mov 2003). The jets in both systems appeared as satellite 
emission components on both sides of the strong HI Balmer 
and Hel emission lines, and they looked similar to the addi- 
tional emission bumps in RS Oph (Fig. 2b) . Similar profiles 
have been observed in the symbiotic nova RX Pup (Mikola- 
jewska et al. 2002) and AR Pav (Quiroga et al. 2002). 

At quiescence, the hot components of RS Oph and 
RX Pup show activity (high states) characterized by the 
appearance of B/A/F shell- type features and variable 'false 
atmosphere' together with the broad and complex emission 
lines and are related to the accretion flow/accretion-driven 
phenomena in these binaries. Other systems that show sim- 
ilar behaviour are CH Cyg and MWC 560 (Mikolajewska et 
al. 2002, 2003). 

Here, we will consider three possible origins for the 
broad component in RS Oph: 

(1) The first possible origin is related to ejection of 
blobs of matter. These blobs could be expelled by a rotat- 
ing white dwarf magnetosphere (Tomov 2003 and references 
therein) or by a jet mechanism. The de-projected velocity 
of the ejection is ~1500-4000 km s _1 , having in mind V em m 
from Table 1 and assuming it is realized normal to the or- 
bital plane, and also the orbital inclination is 30° —40°. This 
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velocity is practically the same as ejection velocities of up to 
3800 km s~ seen during the 1985 nova outburst (Bode 1987 
and references therein). 

(2) Variable disk winds with similar (to RS Oph and 
CH Cyg) velocities were detected in a few cataclysmic vari- 
ables and the whole wind can even turn on and off. The 
terminal velocity of the wind in the BZ Cam system is 
v t ~3000 km s _1 and the time scale of the variability is 30- 
40 min (Ringwald & Naylor 1998). Indeed, their profiles bear 
a noticeable similarity to those in RS Oph at some epochs. 
For Q Cyg the terminal velocity is v t ~1500 km s _1 and the 
events last about 1.5 hours (Kafka et al. 2003). These time 
scales are similar to those of the Balmer line variability of 
CH Cyg but not to RS Oph. It is possible that the disk wind 
of RS Oph varies on longer time scales. 

(3) If the variable broad component of RS Oph orig- 
inates in an asymmetric disk, then a Keplerian velocity of 
1000 km s _1 requires a distance from a 1.4 Mq WD of about 
2xl0 10 cm. At that distance the Keplerian period is 20 min, 
which is considerably different from the observed time scale 
of Ha variability of RS Oph. 

All of the mechanisms considered are related to the loss 
of mass by the system. Whenever such mass-loss occurs, a 
question about its quantitative estimate arises. A quanti- 
tative estimate, however, is possible in those cases where 
the phenomenon can be related to some geometrical model. 
However, we have no geometrical model of loss of mass which 
can be related to the irregular form of the profile of the 
broad component. The velocity distribution of this line pro- 
poses movement of discrete regions (blobs of matter) in the 
emitting environment (Fig. 2). That is why we tried to ob- 
tain a rough estimate of the mass of one "average" blob on 
the basis of its emission and supposing that it is optically 
thin and has also a constant density within some limits. 
From Fig.lb and Fig. 2, the luminous flux in Ha of such 
a blob is about 3xl0~ 13 erg cm' 2 s _1 (W« 3A). This flux 
was corrected for interstellar reddening using E(B— V)=0.73 
(Snijders 1987; Hachisu & Kato 2000, 2001) and the extinc- 
tion law by Seaton (1979). A corrected flux of 1.602xl0~ 12 
erg cm -2 s _1 was obtained. 

For calculation of the emission measure of this blob, 
data about the distance to the RS Oph system, the Ha 
recombination coefficient and the helium abundance are 
needed. For the distance to the system we adopted 0.6 kpc 
according to Hachisu & Kato (2000, 2001). It is not possible 
to obtain the electron temperature and the electron den- 
sity in the region where the Ha broad component is emit- 
ted from observation since we have no indication about the 
appearance in this region of certain lines giving informa- 
tion for those parameters. We will suppose that the electron 
temperature is 20000 K and the electron density is com- 
paratively high being in the range 10 8 - 10 10 cm -3 . Then 
we used recombination coefficients of 5.956xl0~ 14 cm 3 s _1 
and 6.336x 10~ 14 cm 3 s _1 for case B, corresponding to these 
temperature and densities (Storey & Hummer 1995). We 
also adopted a helium abundance of 0.1 which is thought to 
be typical of symbiotic nebulae (Vogel 1993; Vogel & Nuss- 
baumer 1994). 

For calculation of the emission measure we need to know 
the state of ionization of helium in the emitting region. The 
HeIIA4686 line is absent in the spectrum (Dobrzycka et al. 
1996) or is weak (Fig.l) during the quiescent state of the 



system. This means that singly ionized helium is dominant 
in the circumbinary nebula during this state. According to 
the suppositions presented however, the broad Ha compo- 
nent is emitted in a region placed in the close vicinity of the 
hot object in this system where probably helium is mostly 
doubly ionized. That is why here we will assume the state 
of ionization to be He ++ in the emitting region. 

Using a flux of 1.602xl0~ 12 erg cm~ 2 s _1 and den- 
sity of 10 s - 10 10 cm~ 3 , we obtain the emission measure of 
one "average" blob of matter (with a constant density) of 
3.19 x 10 56 - 3.00 x 10 56 cm" 3 . We can also calculate the mass 
of the blob, adopting the parameter fi of 1.4 (Nussbaumer & 
Vogel 1987), determining the mean molecular weight (iraH 
in the nebula. The mass is therefore obtained as 3.8xl0 -9 - 
3.5X10" 11 M Q . [Note that Bode (1987), gives d= 1.6 ± 0.3 
kpc from a variety of measures. If we adopt this longer dis- 
tance, we obtain mass of the blob 2.7xl0~ 8 - 2.5xl0~ 10 
M ]. 



5 CONCLUSIONS 

We report that the Ha emission of RS Oph was strongly 
variable during our 2004 observations on time scales of 
~1 month. No variability was detected on time scales 
of ^1 day. The line was always very wide at its base 
(FWZI?a4600 km s" 1 ) on all spectra from 1986 to 2004, 
consisting of narrow and broad emission components. Vari- 
able emission is detected at velocities up to ±2000 km s _1 . 
Most probably this variable emission is due to either blobs 
ejected from the white dwarf or a variable accretion disk 
wind. The approximate mass of one "average" ejected blob 
of matter, suggested by the variability of the broad emission 
component is 3.8 xl0~ 9 - 3.5 xlO -11 M . We also detected 
variability of HeIIA4686 line on time scales <1 day. To un- 
derstand more fully the nature of the emission line variabil- 
ity of RS Oph we need to acquire a set of spectra with time 
resolution from minutes to days. 
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